Dirac nodal pockets in the antiferromagnetic parent phase of FeAs superconductors 
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We show that previously measured small Fermi surface pockets within the antiferromagnetic 
phase of SrFe2As2 and BaFe2As2 are consistent with a Dirac dispersion modulated by interlayer 
hopping, giving rise to a Dirac point in /c-space and a cusp in the magnetic field angle-dependent 
magnetic quantum oscillation frequencies. These findings support the existence of a nodal spin- 
density wave in these materials, which could play an important role in protecting the metallic state 
against localization effects. The speed of the Dirac fermions in SrFe2As2 and BaFe2As2 is found 
to be 14—20 times slower than in graphene, suggesting that the pnictides provide a laboratory for 
exploring the effects of strongly interacting Dirac fermions. 
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Despite broad similarities in their superconducting 
phase diagrams [1 , a crucial difference lies in the metal- 
lic character of the parent antiferromagnetic phase of 
the iron arsenide (FeAs) superconductors unlike the 
cuprates [2 , 3, 4 . An unresolved question in these mate- 
rials concerns whether the absence of a Mott insulating 
regime is the product of greatly reduced electron corre- 
lations compared to the cuprates [H |5]. An alternate 
possibility is a nodal dispersion that protects a metal- 
lic density of states at the Fermi energy against com- 
plete gapping arising from a strongly correlated spin- 
density wave [6 . Recent angle-resolved photoemission 
spectroscopy (ARPES) measurements [7 reveal a conical 
Dirac-like dispersion within the antiferromagnetic phase 
of BaFe2As2 (see Fig. [l^), presenting a potential route 
for realisation of the latter case. Such a Dirac-like dis- 
persion signals degeneracy in the electronic structure- 
like graphene [8 , suggesting potential topological quan- 
tisation. 

While in two-dimensional (2D) graphene a precise tun- 
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FIG. 1: Portion of the bandstructure of antiferromagnetic 
BaFe2As2 (a) indicating Dirac points (shifted relative to 
/i = 0) responsible for [3 and 7 pockets shown in b. Calcu- 
lation provided courtesy of M. Johannes using the procedure 
detailed in Ref. 



ing of the chemical potential (/i) is required in order to 
access the Dirac node, hopping along the interlayer di- 
rection in the quasi-2D FeAs materials provides an in- 
trinsic tuning to access pockets of fermions positioned 
at the Dirac nodes. Degenerate electron and hole pock- 
ets are expected at the Dirac nodes in such quasi-2D 
systems, evincing characteristic topology. In this pa- 
per, we revisit the small (5 and 7 Fermi surface pockets 
observed via quantum oscillations in the antiferromag- 
netic phase of the FeAs superconductors (SrFe2As2 and 
BaFe2As2) [2l[3] and investigate their origin. In contrast 
to graphene where Dirac points arise within the origi- 
nal band structure, we look for potential association of 
the small Fermi surface pockets with Dirac points occur- 
ring as a consequence of spin-density wave (SDW) band 
reconstruction in the FeAs materials [6 . We probe for 
consistency of the measured jS and 7 pockets (bandstruc- 
ture calculation shown in Fig. [T]) with almost degenerate 
electron and hole pockets located at a Dirac node. By 
a comparison of experimentally measured Fermi surface 
topology and effective mass with predictions, we show 
that the small (5 and 7 pockets are consistent with a Dirac 
dispersion with characteristic speed c* ~ 5— 7 x 10^ ms~^ 
matching that reported in ARPES experiments [7|. 

Figure [2^ depicts how a small interlayer hopping of 
t ~ 10 meV would cause the apex of the Dirac dispersion 
to intersect the Fermi energy even for a small chemical 
potential, giving rise to electron and hole pockets as a 
function of the interlayer dispersion. Vertically stacked 
electron and hole pockets of the form depicted in Fig. [2^ 
touch at their extrema, representing the Dirac point in 
/c-space. By considering a minimal model capturing a 
Dirac electronic dispersion, both the (5 and 7 pockets in 
SrFe2As2 and BaFe2As2 can be simultaneously explained 
by just two parameters (in addition to c*). The un- 
usual lemon-shaped topology of Fermi surface (illustrated 
in Fig. ^p) causes the extremal cross-sectional areas to 
exhibit cusp-like dependences on the orientation of the 
magnetic induction B, accounting for the reported up- 
turn in the magnetic field angle-dependence of the quan- 
tum oscillation frequencies [2*, "3^ . Further predictions are 
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orbits lie in the {kx, k' ) plane, with loci 




FIG. 2: a, Fermi surface schematic according to the disper- 
sion given by Eqn. [l] showing electron and hole pockets with 
the planes of the extremal orbits shown for an arbitrary angle 
between ky and ky. The magnetic induction (not shown) is 
orthogonal to the planes of the orbits, b, The shape of the or- 
bit for several different angles 0, indicating the collapse in the 
radius of curvature /cr at its furthest extremity approaching 
the Dirac point as ^ ^ 90°. 



made for the angle-dependences of the cyclotron effective 
masses and magnetic breakdown tunneling effects, which 
we propose to provide an opportunity for studying exci- 
tations of strongly interacting fermions close to a Dirac 
point. 

We begin by considering a dispersion of the form 



€ = ±hc*\k\ + 2tcos[cA:;,/2] + /i. 



(1) 



where c is the bilayer spacing in the body-centered tetrag- 
onal crystal structure and a non-zero chemical potential 
/i causes asymmetry in size between electron and hole 
pockets [2 , 3 . We neglect the anisotropy in c* identified 
in ARPES measurements [7 and band structure calcu- 
lations [3], so that k = ^k^ + k'^ where k^ and ky lie 

within the Fe As-planes [9 . When B is aligned parallel to 
the crystalline c-axis and /c^, the extremal cross-sectional 
areas and cyclotron effective masses are 



{2t±\^i\f and 
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2t±\^i\ (2) 



and 



refer to 



respectively, where the superscripts 
two different pocket sizes. 

To determine the extremal Fermi surface cross-sections 
for an arbitrary orientation of B, we introduce auxiliary 
coordinates {k'y^ k'^) = Ile{ky^ kz)^ where Yie is a rotation 
matrix and is the angle between B||/c^ and the crys- 
talline c axis. Under such a transformation, the extremal 
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obtained by setting A:^ = and substituting ±/i in place 
of /i in Eqn ([T]). The cross-sectional areas and cyclotron 
effective masses are then obtained by evaluating 



^k,6 



?R:[k^]dky and mf 
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at £ = 0, where 2/cl = 2 cos~^[±|/i|/2t]/c is the length of 
the pocket along kz. 

We first turn to the topology of the /3 and 7 pockets and 
make a comparison with vertically stacked electron and 
hole pockets expected for a Dirac-like dispersion. The ex- 
perimentally measured angular dependence of the /3 and 
7 frequencies in both SrFe2As2 and BaFe2As2 is com- 
pared with predictions of the model in Eqn 1. The ratios 
t/c* and /i/c* in the model can be optimised such that 
the 6>-dependences of the cross-sectional areas A^^ eval- 
uated numerically (lines) using Eqn. [4] in Figs. 3^ and 
b reproduce reasonably well the experimentally observed 
/3 and 7 frequency angle dependences (circles) in both 
SrFe2As2 and BaFe2As2. The extremal cross- sectional 
areas are converted to quantum oscillation frequencies 
by use of the Onsager relation Fq = {fi/2iie)A^Q. The 
ability of the model to reproduce simultaneously the ex- 
perimental angular dependences of two frequencies with 
only two parameters (i.e. t/c* and /i/c*) represents a 
considerable improvement over the prior ellipsoidal ap- 
proximation, which required four parameters in total (i.e. 
two parameters for each frequency [10] ) . The close agree- 
ment with experiment signals the (3 and 7 pockets to be 
consistent with a single cone of Dirac fermions for which 
the Fermi energy intersects the Dirac point at a location 
in k space. 

The next comparison we make is of the measured cy- 
clotron effective masses with model predictions. The 
measured effective mass of the /3-orbit (corresponding to 
A'^ q) in SrFe2As2 is larger than that of the 7-orbit (cor- 
responding to A'^q) (see Fig.pp, [11 ). As anticipated for 
(3 and 7 orbits arising from tne same Dirac cone, the cy- 
clotron effective mass is seen to be ^-dependent (un- 
like for a quadratic dispersion) — (i.e.) higher for the 
larger cross-section (5 orbit. A comparison of the mea- 
sured cyclotron effective masses with simulations shown 
in Fig.jSj: and d can yield individual estimates for all three 
parameters in the model (i.e. t, /i and c* in Eqn. ([T])). 
We use the effective mass data currently available at lim- 
ited angles [2, 3 to obtain these estimates, presented in 
Figs.Js]^ and d. The estimate we obtain for c* (listed in 
Fig. IsT provides further support for the Dirac picture, 
proving consistent with ARPES experiments [7]. 

Given signatures of the validity of the Dirac dis- 
tribution in antiferromagnetic SrFe2As2 and BaFe2As2 
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FIG. 3: Lines indicate simulated field angle-dependences of 
the magnetic quantum oscillation frequencies in SrFe2As2 (a) 
and BaFe2As2 (b) using Eqn. Q and the Onsager relation, 
and the corresponding effective masses (c) and (d). The cor- 
responding simulation parameters (c* = 5.2 x 10^ ms~^, 
t = 9.9 meV, /x = 4.1 meV and c = 6.15 A for SrFe2As2 



and c* 



7.4 X 10^ ms" 



15.3 meV, fi — 5.6 meV and 



c = 6.48 A for BaFe2As2) which determine the form of the 
dispersion given by Eqn. [l] are adjusted to match the exper- 
imental values (represented by circles) from Refs. ^2, |3, .llj 
and Fig. [5b. 



(Fig. [3|, further magnetic quantum oscillation measure- 
ments are indicated to access the quasiparticle excita- 
tions close to the Dirac point. A test of a Dirac disper- 
sion is typically provided by experimental verification of 
the effective mass vanishing as the quasiparticle trajec- 
tory approaches the Dirac point. In graphene, this is 
achieved by using electric fields to depopulate the pocket 
of carriers p!2^. The equivalent depopulation in SrFe2As2 
and BaFe2As2 would require fi to be tuned by doping or 
pressure. 

Angle-dependent measurements provide an alternative 
means for accessing the Dirac point that does not re- 
quire doping. As the angle 6 is increased, the radius of 
curvature /cr (i.e. the radius of a circle tangent to the 
orbit and having the same curvature) at the furthest ex- 
tremity of the orbit (shown schematically in Fig. ^p) un- 
dergoes a dramatic reduction in size, collapsing to zero 
as 6> ^ 90°. The high extremal quasiparticle velocity 
= c*cos^ + (tc/h) sin sm[{cky/ 2) sin 0] turning on a 
small radius translates to a vanishing contribution of a 
quantity we term the 'extremal mass' tue = hk^/vE to 
the cyclotron mass at the farthest extremity of the orbit. 
Although rriE cannot be isolated in quantum oscillation 
experiments made at a single angle, the orbit ally aver- 
aged cyclotron mass given by Eqn. [4] is strongly affected 
by rriE- A consequence of tue vanishing as 6> ^ 90° is 
a marked reduction in the ratio mf/A^ near ~ 90° 
in Fig. |4] The predicted dip in this ratio near = 90° 
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FIG. 4: Simulated ratio (solid lines) of the cyclotron effec- 
tive mass to the frequency (given by the cross-sectional 
area q) for SrFe2As2 (a) and BaFe2As2 (b), yielding a pro- 
nounced dip at ^ = 90° . The dip is a consequence of the van- 
ishing contribution to the cyclotron mass from the extremity 
of the orbit closest to the Dirac point in /c-space. For com- 
parison, we show the same ratio expected for an ellipsoidal 
Fermi surface comprising conventional Landau quasiparticles 
(dashed lines), which is independent of 0. 



(solid lines in Fig. |4| is in contrast to the constant ra- 
tio rriQ/A = mo/Ao expected for an ellipsoidal pocket 
of revolution comprising conventional Landau quasipar- 
ticles (dashed lines). Further experiments accessing the 
cyclotron mass as a function of angle are anticipated to 
provide a confirmation of the experimental distinction 
between massless Dirac quasiparticles in the antiferro- 
magnetic FeAs family and conventional Landau quasi- 
particles. 

Finally, we note that once = 90°, the Fermi sur- 
face topology no longer supports closed extremal or- 
bits. The quasiparticles instead move continuously be- 
tween the electron and hole pockets via the Dirac point. 
Magnetic breakdown tunneling further ensures that open 
orbits will occur for a finite range of angles near 90°. 
The probability = ex.-p{—Bo/B) of magnetic break- 
down (across the gap shown schematically in Fig. |5|i) 
depends on th e charac teristic magnetic breakdown field 
Bq = (7rfi/e) [13 for this Fermi surface topol- 

ogy. This in turn depends on the band gap kc at the 
point where magnetic breakdown breakdown occurs (i.e. 
the separation between electron and hole orbits for a 
given value of 6>, see Fig. [S^i) and on the radius of curva- 
ture of the orbit /cr— both of which depend on the form of 
the dispersion. The dispersion in the immediate vicinity 
of the Dirac point for small kz — ki, has the approximate 
form 6 ^ ±fic*|/c| + hvE{kz — A:l), from which we esti- 
mate kc ^ 2{c* /ve) ki, cot and /cr ~ {ve / c*)ki, cot 
for a point of magnetic breakdown located at a distance 
ky = ki, cot away from the z-axis in Fig. l5k. Substi- 
tuting kQ and /cr into the above expression for Bq^ we 
obtain 



Bo^ { ^]k? cot^ 0. 

e \veJ 



(5) 



Such a strongly angle-dependent magnetic breakdown 
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FIG. 5: a, Schematic of the local geometry in the vicin- 
ity of.the Dirac point where magnetic breakdown occurs as 
90°. b, Amplitude of the 7 frequency in SrFe2As2 
from Ref. after Fourier transforming the raw data hav- 
ing subtracted a third order polynomial, with a fit to 
aT/ smh.(14:.69m^T / B) to obtain the cyclotron mass m^. c 
Amplitude attenuation factor 1 — calculated for the 7 
pocket observed in BaFe2As2 [3]. 

field, which vanishes as ^ ^ 90°, is a unique feature 
of the Fermi surface topology in the vicinity of a Dirac 
point. While open orbits do not contribute to the mag- 
netic quantum oscillations, there will be a loss of ampli- 
tude for the closed orbit given by 1— = 1 — exp(5o/^)- 
The anticipated form of 1 —p^ for the case of the 7 pocket 
of BaFe2 As2 is shown in Fig. [Sj). The dramatic attenua- 
tion of the quantum oscillation amplitude of the 7 pocket 
expected near = 90° in the model may account for its 
reported loss in Ref. [3 at large angle. 

In conclusion, we have shown from the angle- 
dependence of the magnetic quantum oscillations in 
SrFe2As2 and BaFe2As2 that the measured small Fermi 
surface pockets are consistent with Dirac nodal pockets 
arising from a Dirac dispersion with interlayer warping. 



We propose further angle-dependent experiments that 
can be performed to probe the massless nature of the 
quasiparticles approaching the Dirac point, and find that 
the quantum oscillations ultimately become attenuated 
at 6> = 90° due to magnetic breakdown tunneling giving 
rise to open orbits. 

Our finding of Dirac fermions in SrFe2As2 and 
BaFe2As2 implies that the elementary excitations and 
orbital quantization effects are likely to be different from 
those in conventional spin-density wave materials— be- 
coming in fact more like those in graphene ^4j. An 
intruiging aspect of the FeAs materials is not only the 
preservation of a metallic density of states despite a 
large spin density wave gap, but in fact the creation of 
Dirac nodes as a consequence of spin-density wave fold- 
ing, thereby defeating complete Fermi surface gapping. 
The multi-orbital character of the FeAs family appears 
crucial in this creation of Dirac nodes at the meeting 
point of electron and hole dispersions of different orbital 
symmetry [6^. The virtual immunity of the observed 
pockets of Dirac fermions to localization effects could 
potentially explain the persistence of metallic behavior 
despite the very small density of carriers and strong cor- 
relations in the undoped antiferromagnetic materials. Of 
additional special interest is the greatly reduced rela- 
tivistic speed c* ~ 5—7 x 10^ ms~^ in these materials 
compared to 10^ ms~^ in graphene, signalling the sig- 
nificantly stronger correlations of the Dirac fermions in 
SrFe2As2 and BaFe2As2. A pertinent question concerns 
whether strongly interacting Dirac fermions can tunnel 
through barriers or support superconducting pairs over 
long distances (by way of the proximity effect) in a sim- 
ilar manner to those in graphene [14]. 

This work is supported by the US Dept. of Energy, 
the National Science Foundation and the State of Florida 
through the National High Magnetic Field Lab, Trinity 
College (U. of Cambridge), and the Royal Society. We 
acknowledge helpful discussions with M. Johannes. 



[1] S. A. Kivelson & H. Yao, Nature Mater. 7, 927 (2007). 
[2] S. E. Sebastian et al, J. Phys.: Cond. Matt. 20, 422203 
(2008). 

[3] J. G. Analytis et a/., Phys. Rev. B 80, 064507 (2009). 

[4] G. F. Chen et aL, Chinese Phys. Lett. 25, 3403 (2008); 
M. Rotter, M. Tegel, & D. Johrendt, Phys. Rev. Lett. 
101, 107006 (2008); D. J. Singh, Phys. Rev. B 78, 094511 
(2008); J. Zhao, Nature Mater. 7, 953 (2008); J. H. Chu 
et al, Phys. Rev. B 79, 014506 (2009). 

[5] Q. Si et al, New J. Phys. 11, 045001 (2009). 

[6] Y. Ran et a/., Phys. Rev. B 79, 014505 (2009). 

[7] P. Richard et al, < lhttp://arxiv.org /abs/0909 .0574[ >. 

[8] M. L Katsnelson, Mater. Today 10, 20 (2007). 

[9] Anisotropy in c* can split the observed f3 and 7 frequen- 
cies into two sets of frequencies for an arbitrary orienta- 
tion of B, which may be too small to be resolved experi- 



mentally in high magnetic field studies. 

[10] In the ellipsoidal Fermi surface model [2j [3], each fre- 
quency requires two parameters— the cross-sectional area 
at ^ = and the ellipticity, totaling 4 for both pockets. 

[11] The f3 pocket in BaFe2As2 [3] is reported to have an 
anomalously low mass and large error bar, perhaps re- 
flecting its small amplitude in Ref. [3 . 

[12] K. S. Novoselov et a/., Nature 438, 197 (2005). 

[13] D. Shoenberg, Magnetic oscillations in metals (Cam- 
bridge University Press, Cambridge, 1984). 

[14] A. H. Castro-Neto et aL, Rev. Mod. Phys. 81, 109 (2009). 



